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Piercing the Vainshtein screen with anomalous gravitational wave speed: 
Constraints on modified gravity from binary pulsars 
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By using observations of the Hulse-Taylor pulsar we constrain the gravitational wave (GW) speed 
to the level of 10“^. We apply this result to scalar-tensor theories that generalize Galileon 4 and 
5 models, which display anomalous propagation speed and coupling to matter for GWs. We argue 
that this effect survives conventional screening due to the persistence of a scalar field gradient inside 
virialized overdensities, which effectively “pierces” the Vainshtein screening. In specific branches 
of solutions, our result allows to directly constrain the cosmological couplings in the effective field 
theory of dark energy formalism. 


Introduction - Modifications of General Relativity 
(GR) that explain the acceleration of the Universe can 
display a gravitational wave (GW) speed ct 1 (we use 
units h = c = 1). What are the observational constraints 
on this parameter? In some given model, ct can be ex¬ 
pressed as a specific function of the (post-Newtonian) 
parameters of the theory, and thus constrained indi¬ 
rectly with solar system tests (see e.g. i)- On the other 
handjcosmological observations limit ct to the 10% level 
{e.g. IH). In Ref. 0, Moore and Nelson observe that sub¬ 
luminal GWs would be Cherenkov-radiated by particles 
traveling faster than ct. By looking at high energy cos¬ 
mic rays data, the authors manage to constrain this effect 
to the impressive level of 10“^^. We notice, however, that 
the typical energy of the corresponding radiated gravi¬ 
tons, 10^° GeV, is well above any reasonable cut-off of 
the modified gravity theories for cosmic acceleration. It 
is not difficult to envision e.g. Goldstone modes in spon¬ 
taneous Lorentz breaking situations that are subluminal 
at low frequencies and recover relativistic propagation 
above the symmetry breaking scale . With binary pul¬ 
sars timing data, in this letter we obtain for ct looser 
limits (^ 10~^), which however apply to frequencies that 
are relevant for an effective theory of dark energy. 

One obvious objection is that scalar tensor theories 
generally come equipped with screening mechanisms, al¬ 
lowing to recover the stringent tests of gravity in the 
galaxy and in the solar system. Among these, the Vain¬ 
shtein screening Q is particularly efficient, and relevant 
for those scalar tensor theories that display anomalous 
GWs speed. What screening guarantees, however, is 
the suppression of the contribution of the scalar field cj) 
to the total gravitational attraction between bodies in 
the Newtonian approximation. In a screened situation, 
the fluctuations of the metric field—gravitons—are left 
as the only mediators of long-range interactions. But 
not necessarily do they behave as in GR. The point is 
that the background value of the scalar (po , although not 
directly participating in gravitational interactions, gen¬ 
erally maintains a non vanishing gradient that sponta¬ 
neously breaks Lorentz symmetry. In such a situation, 
the effective gravitational Lagrangian need not be that 
of GR, even if it involves only massless gravitons. The 
Vainshtein screen is •pierced. 


The same mechanism is responsible for other violations 
of the screening considered in the literature. In simple 
cases, deviations from GR boil down to a spacetime vari¬ 
ation of the Newton constant Gn- Refs. @,0 use Lunar- 
laser-ranging to constrain this effect, obtaining limits on 
modified gravity models that are comparable in size to 
those obtained here. Preferred-frame effects 0 and pos¬ 
sibly anomalous values of the gravitational slip parame¬ 
ter 7ppn 0 , [13 (see also the following on this) have also 
been discussed in the literature. 

Quadrupole formula, revisited - For the sake of 
generality, we will consider a two-fold modification of GR 
encoded in the following Lagrangian for the GWs sector: 


£ = 


1 

GJTrLrg,^ 


E 




( 1 ) 


where -I-, x represent the two polarizations of the GWs. 

First, we allow for a coupling of GWs to matter, Ggw, 
possibly different than the Newton’s constant Gm in¬ 
ferred in the Newtonian limit via the Poisson equation. 
Indeed, in addition to the radiating gravitons described 
by the above Lagrangian, we have the potential gravi¬ 
tons [ll| , responsible for the bound of the binary system. 
In modified gravity theories with an additional scalar de¬ 
gree of freedom, the scalar sector also becomes radiative. 
We will rely on Vainshtein screening while assuming that 
the contribution of the radiated scalar to the variation of 
the binary system period is negligible, as was shown to 
be the case in specific models [l^ . 

The second modification that we consider is that GWs 
can propagate at a speed cp different from the speed of 
light. We assume here that such a speed is constant, 
direction- and polarization-independent. This statement 
is exact in the limit of a constant gradient for the back¬ 
ground scalar field cpQ, and in the reference frame where 
such a gradient is along the time direction. In the follow¬ 
ing we will quantify the corrections due to the presence 
of a spatial component of the gradient, and argue that in 
realistic situations such a component is negligible. 

It is interesting to revisit, step by step, the standard 
derivation of the quadrupole formula {e.g. 0 ) at the 
light of these modifications. First, we want to estimate 
the energy flux of a GW across a spherical surface at large 
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distance r from the source. The standard expression can 
be modified, essentially, by dimensional analysis {e.g. by 
rescaling the time as dt = crdt' so that the GR formulae 
can be applied straightforwardly). We find 

d F f 

Id = 327rcTGgw J ^ 

where (...) means average over a region of spacetime 
much larger than the GW wavelength. On the other 
hand, at the lowest (quadrupole) order in the velocity 
expansion, the radiated amplitude of GWs from a given 
source is obtained with the usual formulae, barring the 
replacement Gn —^ Ggw and the different retarded time 
at which the source is evaluated, 

[lij]quad = ^ (3) 


where is the transverse-traceless projection of the 
quadrupole moment Qij of the source. Note that 
appears in © after using the energy-momentum conser¬ 
vation of matter. Since the matter sector has the usual 
Lorentz symmetry, the time derivatives acting on Qjj^ 
do not introduce any additional factors of ct- As in 
the standard calculation, the way such a projection is 
made depends on the direction of the GW and this should 
be taken into account when calculating the surface inte¬ 
gral ©. This results in the following total power emitted 

Pquad=^(Q.A,). (4) 


This expression coincides with the formula obtained in 
0 for Horava gravity. 

Binary pulsar constraints - By the above modi¬ 
fied quadrupole formula, binary pulsars observations will 
allow us to constrain the combination crCgw, modulus 
some assumptions on the expressions of the Keplerian 
parameters of the bound system that we detail in the 
following. The emission of GWs results in a decrease of 
the orbital period Pt [l^. Mutatis mutandis, we get 


Pb = - 
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mp-mcinip + me) 


where e is the eccentricity of the Keplerian orbit and 
rup and me are the masses of the pulsar and its compan¬ 
ion, and we have temporarily reintroduced (just here and 
in (O) the dimensional speed of light c. As explained be¬ 
low Eq. o, we assume potential gravitons and radiative 
gravitons to couple to matter with different strengths. 
Note the different roles in the derivation played by Gn, 
coming from the formula of the orbits, and Ggw, coming 
from the actual emission of gravitational waves. 

We use the most accurate available data on P{,, those 
of the Hulse-Taylor pulsar (PSR B1913-I-16) [l5|, with 


Parameter 

Description 

Value 

e 

eccentricity 

0.6171334(5) 

P 6 (days) 

period 

0.322997448911(4) 

w(deg/yr) 

periastron advance 

4.226598(5) 

7 (ms) 

Einstein delay 

4.2992(8) 

Pb 

period decay 

-2.423(1) X 10"^^ 


TABLE I: Orbital parameters for PSR B1913-I-16 from [l^ . 



FIG. 1: Mass-mass diagram for PSR B1913+16 (the Hulse- 
Taylor pulsar) based on the post-Keplerian parameters w 
(black), 7 (red) and Pb (blue). Varying the combination 
ctG^w/Gn amounts to shifting the 1 -cr stripe of Pb. 


the orbital parameters shown in Table . Before using 
this information, we need the standard expressions for 
the advance of the periastron ui and the amplitude of the 
Einstein delayy [T^ , which also depend on the Keplerian 
parameters e and Pb, and on the masses Wp and me- 
We can thus use the binary pulsar data to constrain the 
combination c^Ggw, in addition to the two masses. 

While the expression of 7 is derived, essentially, in the 
Newtonian approximation, a comment regarding the pa¬ 
rameter oj is in order here. In a modified gravity set-up, 
such a quantity depends on both post-Newtonian param¬ 
eters 7 ppn and ,5ppN [ 110 . Since we are aiming (see 
below) to a precision of 10 “^, we rely on solar system 
tests, which constrain yppN and /3ppN at the levels of 
10 “® and 10 “® respectively, and uj directly and indepen¬ 
dently, for Mercury, at the level of I0“® (isl. 1^. 0 . 

We can now proceed to construct the mass-mass di- 


lEq. ® is calculated in the orbiting system reference frame 
which is accelerated with respect to the solar system barycenter 
frame [la- This effect, known as Shklovskii effect, gives an extra 
APfj^gai = —0.027 zb 0.005 x 10“^^ which should be subtracted. 
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agram and the corresponding constraints on ctG^^ as 
shown in Fig. [T] We see that the binary pulsar data 
meet in a small region of the (m^, rric) plane. GR predic¬ 
tions fall in the intersection of uj and 7 within about Icr 
confidence. By imposing compatibility of the three con¬ 
straints at the ItT level we obtain the following bound; 

0.995 < < 1 . 00 . ( 6 ) 

Gn ct 


Symmetries and scalar field gradients - We are 

now to discuss the implications of the above bound on 
concrete scalar tensor models for dark energy. First, it is 
helpful to consider the basic structure of the simplest sce¬ 
nario that displays an anomalous GW speed, the quartic 
galileon model, with Lagrangian 


Cgal _ 
i-4 — 


'as 


(□^)2 




(7) 


In the above, X = df^cfid^cj) and A is some energy scale 
of the order of A ~ with Hq the Hubble pa¬ 

rameter today. By inspection of the second term inside 
the square brackets, we see that the covariant derivatives 
generate a term quadratic in the Christoffel symbols. In 
the presence of a background field (f>o with non-vanishing 
timelike gradient, such term contributes to the quadratic 
Lagrangian for the gravitons hij as ~ , thus modifying 

the propagation speed of GW. It is immediate to see that 
Ct is dependent only on the gradient of (j)o in this case. 
Only when Xp(j)o vanishes does ct go to one. 

But for theories enjoying shift symmetry </>—>■(/> -|- 
const., the actual value of the scalar is irrelevant, and 
there is no evident mechanism for it to detach from cos¬ 
mic evolution and become constant inside a virialized ob¬ 
ject. This means that we do expect, in general, a non¬ 
vanishing scalar gradient inside screened environments— 
a “local remnant” of the expansion of the Universe. 

Indeed, the profile of a cosmologically evolving scalar 
field in the presence of a matter source is easily estimated 
for those theories that enjoy a further, galileon, symme¬ 
try which makes a constant gradient of </>, and not 
only its actual value, irrelevant. In Minkowski space, this 
is defined as the invariance under Xpfi —>■ -|- bp, with 

bp a constant vector. Let (j>Q°^’^^{t) be the cosmological 
solution obtained under the assumption of homogeneity 
and isotropy. Well inside the Hubble radius, where the 
metric is similar to Minkowski, this is effectively a field 
configuration of constant gradient. Once we find a suit¬ 
able radial solution (r) vanishing at infinity around 
some localized matter source, in virtue of galileon sym¬ 
metry, we can simply add the two solutions. 


Mr,t) c, + ( 8 ) 


Galileon theories are a combination of 5 Lagrangian 
terms with an increasing number of fields (j). 

As a case study, let us consider the quartic galileon (O . 
The cosmological gradient for this theory is given by 


^cosm ^ HoMp ~ A^iL|j"^. On the other hand, the analy¬ 
sis of the Vainshtein mechanism near a spherically sym¬ 
metric object shows that the radial gradient of the scalar 
inside the screened region for quartic Galileon is con¬ 
stant, ^ {M/MpY^^A^ ^ TyA^, where we have 

introduced the Vainshtein radius ry — {M/MpY^^A~^ 
and M is the mass of the matter source. In summary. 


^0 ^ ’ 


(9) 


which shows that a localised source contributes a very 
mild radial component to the total gradient of the field. 
For example, the Sun has ry 1 kpc so this ratio is 
of order ~ I0“®. In comparison, our peculiar velocity 
with respect to the GMB gives a much larger (effect'^ 
10“^). Our estimates are in agreement with the explicit 
numerical calculations of [^ . 

Gravity inevitably breaks the symmetry Vpcj) —> V^^-|- 
bp, if anything, because there is no such thing as a con¬ 
stant vector bp in a general spacetime. However, we can 
apply the above estimates to all scalar tensor theories 
that reduce to galileon in the decoupling limit, formally 
defined as Mp 00 while keeping A constant. Am ong 
these, theories with weakly broken galileon symmetry |24l | 
have their Lagrangians protected against quantum cor¬ 
rections. 

Cosmological EFT operators - We have just shown 
that 4>'q <C ^0 (even) inside the Vainshtein radius, where 
the non linearities in the scalar can become important 
but the metric is very close to Minkowski. The most 
general quadratic Lagrangian for the metric fluctuations 
in the presence of a background scalar field of constant 
timelike gradient is conveniently studied within the effec¬ 
tive field theoryfEFT) formalism for cosmological per¬ 
turbations [25l - [^ . By choosing the time coordinate to 
be proportional to the scalar field (unitary gauge), all 
degrees of freedom are transferred to the metric, cho¬ 
sen to be the one minimally coupled to matter {Jordan 
frame). A limited number of operators capture the linear 
dynamics of the most general scalar-tensor theory with 
an equation of motion of at most second order for the 
propagating scalar fluctuation (^ . Among such opera¬ 
tors, only three affect the pure graviton sector. 


^ M2 


R +€4 {SK^^SKi. 


SK^) - ii ^^'>RSN 


( 10 ) 

where R is the Ricci scalar, is the perturbation of 

the extrinsic curvature of the t = const, hypersur¬ 
faces, ^^^R their Ricci scalar and 6 N the perturbation of 
the lapse function. M, €4 and €4 are time dependent co¬ 
efficients. In GR, M = const., €4 = €4 = 0. The above 
operators arise e.g. in the class of models introduced 
in as a generalization of Horndeski theory, which is 
the most general scalar-tensor theory with equations of 
motion of at most second order [ 2 ^, [ 3 ^ . We refer the 
reader to [3l| for the expressions of €4 and q as func¬ 
tions of the full Beyond Horndeski Lagrangians. 
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To study the effects of the terms (flO)) it is convenient to 
switch to Newtonian gauge on a Minkowski background. 
By forcing a time-diffeomorfism t —t + tt, the fluctua¬ 
tions of the scalar field tt reappear in the action, after 
which we can fix the metric to have the form 


ds^ = —(1 -I- 2<f’)dt^ -I- [(1 — 24')(5y -|- "fij] dxMa;-’ , (11) 


where $ and 4' are the two Newtonian potentials and "fij 
represents the transverse traceless graviton. At highest 
order in derivatives the quadratic Lagrangian reads [3lj | 






^ {ifj - 4(V7*,)") (12) 


T 

- _ 2c-2(i + aff)V$V4'-k 

ciTT^ — C2(V7r)^ -I- mixing terms 


where we have defined the GWs speed = (l-|-e 4 )“^ and 
the beyond Horndeski parameter an = e 4 —£ 4 . In the Jor¬ 
dan frame there is no direct coupling of tt to the matter 
fields, but the scalar-metric mixing terms schematically 
indicated in (HU, of the type V 4 >V 7 r, V'J'Vtt and ’^'k. 
When a n ^ 0 the higher derivative term VIE'Vtt also 
appears [^. The explicit form of the last line of dni) 
depend on all the operators of the EFT— i.e. also on 
those omitted in (flUl) — as well as on the time derivatives 
of £4 and £4, and is responsible for the rich linear phe¬ 
nomenology of dark energy, in which the tt fluctuations 
play a dominant role [32 h34| |. 

However, in the vicinity of a localised matter source the 
TT fluctuations become irrelevant because of the screening, 
and we can thus forget about the third line of (fT^ . As 
long as the on-shell gravitons 7 ^, and the Newtonian po¬ 
tentials $ and 41 can be considered as short-wavelength 
fluctuations on top of a constant background scalar field 
gradient, the first two lines of (fT^ can be borrowed from 
cosmology and applied to general set-ups. This is the 
case for GWs of wavelengths much shorter than the dis¬ 
tance from the source. From the first line of m we 
can read off Ggw = For a given (shift- 

symmetric) theory, the cosmological value of ct (equiv¬ 
alently, of the EFT parameter £4) can be calculated as 
a function of X = —(/)g [ 2 ^, |3l|. If such a gradient ac¬ 
quires a spatial component (f )'—either along the radius 
from a matter source, or in the direction of our motion 
w.r.t. the GMB frame —ct simply transforms as a ve¬ 
locity under a boost of speed v = (j)' /4’ ^-nd becomes 
direction-dependent. Along the two principal directions 
the boosted velocity reads 


astro 

C'T' 


Ct{X) i V 
1 i Ct(X)v 


(13) 


We are left with the second line of which can 

be used to describe the dynamics of the scalar potential 
gravitons in the Newtonian approximation. However, its 
applicability to general screened situations is more sub¬ 
tle. Since the Newtonian potentials and the background 


field (j)o are generated by the same source, they are of the 
same typical wavelengths, and the constant gradient ap¬ 
proximation for (jjQ is not guaranteed to work. By substi¬ 
tuting ~ —24)pm, one would obtain the relation 

between the two Newtonian potentials 7 ppn = 4t/4) and 
the Newton constant by the Poisson equation: 


7PPN = 


1 + Off 


Gn = 


87rM2(l -I- auY 


(14) 


The study of spherically symmetric configurations in 
the full beyond Horndeski models confirms that the 
above always correspond to one available branch of solu¬ 
tions d,[l3- Theories with terms up to {type-4) 

show a total of three branches, in two of which the GR 
result 7 ppn — 1 is recovered inside the Vainshtein ra¬ 
dius. However, for beyond Horndeski of type-5 (terms 
up to terms up to (V^^)^), there appears to be now way 
to recover the GR value in any of available branches. We 
would like to emphasize that the branch corresponding 
to (dH), always present, also develops non-linearities in¬ 
side the Vainshtein radius. A closer inspection of the 
solutions in Q shows, however, that for this branch the 
relevant non-linearities are in the mixed 7r-4> and 7r-4' 
sectors, and not in the self-interactions of the scalar, as 
it is usually assumed. 

The branch recovering yppN — 1, when available, is 
often taken as the appropriate solution inside virialized 
objects [ 3 ^ [ 3 ^, also because it matches the asymptoti¬ 
cally flat solutions in some specific cases [s^. However, 
which branch applies to realistic scenarios is ultimately 
selected by the time evolution. The point is to under¬ 
stand, case by case, which solution continuously evolves 
from the (unique) linear configuration describing a tiny 
perturbation in the early Universe, and this will depend, 
in general, on the details of the theory. 

Observational constraints and discussion - 
Within the “linear branch” of solutions o, the cos¬ 
mological EFT parameters £4 and an are tightly con¬ 
strained. First, the bound (0) turns into a constraint 
for the combination of parameters (1 -I- an)'^/ ct- At the 
same time, as already noted e.g. in the value of 
the slip parameter in the linear branch (I14|) is power fully 
constrained by the Cassini spacecraft experiment 
7PPN — 1 = (2.1 ± 2.3) X 10“®. This combines with our 
binary pulsar result as in Fig. [2] 

We would like to stress, however, that beyond the de¬ 
tails related to the specific branch of solutions, the bound 
on the GW speed is very general: Hulse-Taylor Pulsar’s 
observations constrains ct at the level of 10“^, barring 
remarkable and unlikely cancellations with the (linear 
and non-linear) physics that determines the orbits of the 
bound system. Our result applies to all dark energy 
models in which gravity is modified enough to display 
a different speed for GWs. Within scalar-tensor theories, 
in particular, we have considered galileon 4 and 5 type 
models and its generalizations and argued that the effect 
is not screened in general, because it is related to the 
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FIG. 2: Combined constraints in the {c^,aH) parameter 
space for the linear branch of solutions (1141) . The (tight) 
bound the Cassini measurement (black curve) and the 
light-blue stripe corresponding to the Hulse-Taylor pulsar 
bound obtained from the top panel. Within Horndeski theo¬ 
ries, because Gn = c^Gg„, the bound turns into a slight 
preference for superluminal propagation. 


the persistence of the (cosmological) scalar field gradient 
even inside conventionally Vainshtein-screened regions. 
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